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Abstract 

We expose briefly the role of entropy in the early universe and in 
particular the importance of searching for new mechanisms of entropy 
production. We describe a mechanism that shows how entropy is 
produced during early annihilations and under which conditions the 
production is not negligible. A MeV r neutrino provides an interesting 
application. 

This paper has been presented in the Nato School Session on 
"Generation of Cosmological Large-Scale Structure" (Erice, Novem- 
ber 1996) and it will appear in the proceedings of the school (ed. 
D.N.Schramm). 



1. Introduction 

Entropy is not a fundamental quantity as e.g. energy density, but it pro- 
vides an auxiliary function useful to describe the physical processes in the 
early universe. This is evident when thermodynamical equilibrium applies 
(quasi-static expansion) because entropy is conserved and thus it represents 
a constant of motion. This is a traditional result || that can be shown 
easily by applying to the cosmological expansion the usual thermodynam- 
ical relation d(pR 3 )/dt = TdS/dt — pdR 3 /dt + J2i PidNi/dt and imposing 
the energy-momentum tensor conservation equation d(pR 3 )/dt = —pdR 3 /dt. 
Then necessarily: 

dS 1 ^ dNi 

But as we are assuming the equilibrium also in the chemical reactions, the 
right sum is zero and hence dS/dt = 0. Entropy conservation allows for 
the straightforward derivation of the R — T relation. If all species are in 
radiative equilibrium then S oc g^iTR) 3 , where Qr = const are the effective 
degrees of freedom, and simply TR = const, otherwise, the definition of 
is generalized to the case in which some species are not ultrarelativistic by 
introducing a function gs{T) that can be easily calculated || and in this case 
RT oc [(?s(T)] _1//3 . The RT factor is important because the particle numbers 
of any species in radiative equilibrium, and in particular the photon num- 
ber iV 7 , is oc (RT) 3 . Thus entropy conservation provides the way to know 
how the photon number changed during the early universe history. During 
the annihilations of some particle species the degrees of freedom decrease 
and photons are produced because the disappearing species releases its en- 
tropy to the radiative plasma including photons. Sometimes in the literature 
this is called entropy production while it is only entropy exchange. Photon 
production is an important quantity to be known because it is needed to 
calculate the relic abundances (n/n 7 ) of particle species or of the barionic 
number. In fact if at some time tf a particle or charge number Nj is frozen, 
we are often able to calculate the abundance (n/n 7 )/ at that time and then 
the relative relic abundance can be expressed as (n/n 7 ) = (n/n y )f(N^/N^). 
This expression shows how photon production dilutes the abundance from 
freezing until the present. For this reason it is convenient to define a dilution 
factor f = N^/N*. Now it is easy to understand the effect of an entropy 
production: it gives a further contribution to photon production and hence 
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to /: / = f g ■ fs (where f g = g } s /g S o and f s = S /S f ). We point out that 
the more recent the entropy production is the more effective it is, because it 
dilutes all densities previously frozen. 



2. Entropy production calculation through the L.W. equation 

Let us consider the classical problem of the freezing of particle abundances 
during annihilations. Some particle species h can annihilate eventually into 
photons through a lighter particle species I: h + h — > I + I — > . . . — ► 
7's. Assuming that the elastic reactions rates are strong enough to support 
the kinetic equilibrium (thus a temperature is defined for the h and it is 
kept equal to photon temperature), the h distribution function will assume 
an equilibrium form: fh(p,t) — f(p,t) = {e^^^^p"^^ + l}" 1 where fih is 
commonly called the pseudo-chemical potential to underline that it does not 
obey, in general, the chemical equilibrium condition Hh + f-h = 0. In this way 
the entropy production rate can be obtained easily from thermodynamics 
using (0), and in our specific case becomes: 

dt T alt U 

This expression can be obtained from statistical mechanics as well. It can be 
noticed that in order to have dS/ dt ^ two conditions must be satisfied: the 
first one is departure from chemical equilibrium (fih + p,^ ^ 0) and the second 
is a decreasing particle number (dNh/dt ^ 0) which means the presence of 
annihilation processes not balanced by pair production processes. If we want 
to perform a quantitative calculation of the entropy production, we must use 
some model to calculate fih and dNh/dt. The simplest way to describe the 
annihilations is represented by the special Lee- Weinberg equation: 



— - A 

dy 



N 2 (y)-Nl(y) (3) 



(y = T/m, N = N h /N in , N in = N h (T > m h ), A = 0.055(g h / ^gf)a mm PL ). 
In the modern literature this equation is usually written using the variable 
Y = n/s (s is the entropy density) where Y oc iV if entropy is conserved. 
But as we are relaxing entropy conservation assumption we must use a vari- 
able like N. In this equation particle-antiparticle symmetry is assumed. 
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This is not restrictive because it means only that the departure from equi- 
librium and the consequent freezing must happen for temperatures when 
(N h — Nj^/Nh <C 1. Otherwise it is well known || that departure from equi- 
librium would be greatly delayed and therefore entropy production would 
be negligible. Thus we can set fih = fih = p> and the chemical equilibrium 
condition becomes fx = 0. Moreover, it is an implicit approximation of the 
equation that z = fx/m — yln[N(y)/N eq (y)]. In Fig. 1 we show the solutions 
N of Lee- Weinberg equation for different values of A and the correspondent 
z . 




Finally we can write the entropy production rate through N obtaining: 

1 dS(y) _ 135£(3) z(y)dN(y) 
St dy 4vr 4 y dy U 

where = ghSm/ 1 9 % s- By integration one can obtain the entropy production 
AS(y,yi n ). The results are shown in Fig. 2. 

The rates present a maximum because entropy production is zero before 
the system goes out of equilibrium (i.e. fx — 0) and after it tends again to 
zero because the annihilation rate decreases (dN/dt — ► 0, while fx reaches its 
asymptotical value m). Moreover, one can notice the existence of a value for 
Aq for which entropy production is maximum. In fact both for Aq — > and 
for Aq — > oo entropy production also tends to zero: in the first case because 
particles would not annihilate at all, though the system is out of equilibrium, 
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and in the second because the system never goes out of equilibrium. Thus 
there must exist a maximum. It is reached when Aq = 3.65 and Xf = 2.5, 
where we defined Xf as that value for which N —N eq = 0.5N eq . The maximum 
production of entropy occurs for a freezing during which /i-particles are semi- 
relativistic. 




Figure 2. 



In the asymptotical limit (x = m/T — > oo) all expressions become simpler: 
dN(x)/dx — > —AqNq/x 2 , z(x) — > 1 and dS(x)/dx — >oc AqNq/x: hence 
asymptotically the entropy production AS" = / dx(dS/dx) does not stop but 
it increases logarithmically, thanks to the relic annihilation. 

3. Estimation of maximum entropy dilution factor 

In order to estimate the maximum dilution that the mechanism can pro- 
duce we must integrate the entropy production over the entire early universe 
history, until matter-radiation decoupling. Another consideration to be made 
is that during the early universe the degrees of freedom decrease, and this 
requires the generalization of the calculation of dN/dt done with the sim- 
ple Lee- Weinberg equation which assumes degrees of freedom and entropy 
constant. At the same time we can also include the correction due to the 
variation of entropy. It is possible to introduce these modifications only in 
the asymptotical regime (which gives the major contribution to the entropy 
production). Hence during this phase the decoupled equations (|3|) and (§) 
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can be replaced by two coupled equations: 



dS 135 £(3) g h ldN dN . g Si * - 2 



A ■ ^^N{ (5) 



dy 4 7T 4 gfydy dy y/gi*S(y) 
(S = S/Sin, g S i = gsi/g l s, 9 P i = 9 pi/ 9™) that yield the solution 



Sdec — S\ — S* 



-,, n a 9h ^ 9Si, Xi-i 135 f (3) - 2 



If the second term in the square root is much smaller than 1 then, ex- 
panding at first order, the logarithmic behaviour is restored. Until now we 
supposed that all components are coupled to photons and therefore fs = Sdec- 
If there is a decoupled component it must be considered that some fraction 
of the entropy can be given to the decoupled component and thus does not 
contribute to the dilution. 

Now we must find the maximum value that the expression (|6|) can assume. 
It depends on two parameters: A and T in (gh can be set equal 2). We have 
already indicated that Aq = 3.65 yields the maximum value for AqNq . To find 
the appropriate T in there are two opposite considerations: to get the greatest 
temperature range requires annihilations as early as possible, but to have the 
greatest fractional weight of /i-particles (= gij 9 % s) requires annihilations that 
happen when as few particle species as possible have survived. The second 
consideration suggests to take a time interval such that g™ is minimum, and 
the first suggests the insertion of annihilations as early as possible inside this 
interval. Thus the best situation is when freezing starts soon after electron- 
positron annihilations so that m ^ m e /2 ~ 0.25MeV and from (^|) fs = 
1.43 is obtained. The second possibility is that freezing starts soon after 
muon annihilations (rrih ^ ~ 50MeV) and in this case fs = 1.31 is 

obtained. 



4. Application: Mev r-neutrino 

In this section we want to consider a realistic case that produces a dilution 
factor as much as possible close to the order of magnitude found earlier. This 
time we must substitute <jq = const with the thermally averaged cross section 
(cannVMc^iy), that can be calculated with a single- integral formula valid 
for any value of Xf. Unfortunately it is not realistic an hundred Kev particle 
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that undergo the semirelativistic freezing (x/ ~ 2) we need, though further 
investigations are in progress considering Majoron or Axino. An Mev r- 
neutrino is a particle that would freeze with the right Xf and it is light enough 
to have a good fractional weight. This then seems to be the best application 
of the mechanism we described and also because of the possible consequences 
that it could have on its mass constraints. In last six years constraints from 
nucleosynthesis || [7|, SN1987A ||, have been greatly improved, reducing 
the possibilities for an MeV neutrino from an astrophysical point of view. 
At the same time the laboratory experiment lower limit has been lowered 
to 24 MeV and thus in future years there could be an important test for 
astroparticle physics. This is why we think that it is worthwhile to continue 
investigations on MeV r neutrino. Here, however, we want mainly to provide 
an application for the general mechanism we described. 
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Figure 3. Entropy dilution factor: MeV Dirac v T (left); MeV Giudice v T (right). 

We performed the calculations for a Dirac neutrino and for a neutrino 
with a magnetic moment <^ 2 • 10~ 9 (m u /10MeV)fiB [10| - In the first case 



it must be considered that v T can annihilate into e + e~, v^v^ v e v e and that 
for T ^ lMeV the neutrinos are decoupled: this reduces of a factor ~ 1/3 
the dilution factor from the maximum and the result is fs = 1.01 4- 1.08 for 
masses m = 2 4- 2AMeV (see Fig. 3 left). In the second case the annihilation 
through photon exchange in e + e~ is dominant. In this case the factor 1/3 is 
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avoided and the dilution factor can approach the maximum value permitted 
(fs = 1.1 -r- 1.3 for masses m — 2 4- 2AM eV as shown in Fig. 3 right). 

5. Conclusions 

In the introduction we clarified the difference between the simple entropy 
transfer from the annihilating particle species to the radiative plasma and 
an effective entropy production. We presented here the possibility to have 
a significant entropy production in the late phase of the early universe if 
the annihilating particle species undergoes a semirelativistic freezing. This 
mechanism can dilute all relic abundances in the early universe by a factor 
1.5 as maximum value. In the last section we presented an application that 
realizes this possibility: an MeV v T . In the general cases of nonrelativis- 
tic and ultrarelativistic freezing, entropy production is a negligible effect as 
usually assumed. 

This talk is based on the work [|TJ . P. Di Bari has been fortunate to have 
V.Berezinsky as the supervisor in his "Tesi di Laurea". He wishes to thank 
A.D. Dolgov for valuable suggestions and M. Lusignoli for his helpful role 
during the work of thesis. He also thanks all people and the director of Gran 
Sasso Laboratories P. Monacelli for their warm hospitality during the last 
two years and Lori Gray for her help with the English. 
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